function [b1DiffvsGS,cumB1DiffvsGS] = diffB1SetAgainstGoldStandard(gvolmask, vol)
%diffB1SetAgainstGoldStandard Calculates the difference of a set of volumes
%relative to a gold standard volume. 
%
%   Detailed explanation goes here

    for ii=1:6
        for jj=1:6
            for kk=1:length(vol)
                result{ii,jj,kk} = sum(sum(abs(gvolmask-vol{kk}{ii,jj})));
            end
        end
    end
    b1DiffvsGS=cell2mat(result);
    cumB1DiffvsGS=sum(b1DiffvsGS,3);
    
end

